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' Abstract. A sequence of rational functions in a variable q is g-holonomic if it satisfies a linear recursion 

JH , with coefficients polynomials in q and q" . We prove that the degree of a g-holonomic sequence is eventually 

a quadratic quasi-polynomial, and that the leading term satisfies a linear recursion relation with constant 
coefficients. Our proof uses differential Galois theory (adapting proofs regarding holonomic D-modules to 
the case of q-holonomic D-modules) combined with the Lech-Mahler-Skolem theorem from number theory. 

Ch , En route, we use the Newton polygon of a linear g-difference equation, and introduce the notion of regular- 

singular q-difference equation and a WKB basis of solutions of a linear g-difference equation at 5 = 0. We 

^^ then use the Skolem-Mahler-Lech theorem to study the vanishing of their leading term. Unlike the case 

^i_i-/ ' of g = 1, there are no analytic problems regarding convergence of the WKB solutions. Our proofs are 

r) , constructive, and they are illustrated by an explicit example. 
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1. Introduction 

1.1. History, q-holonomic sequences appear in abundance in Enumerative Combinatorics; [PWZ96, Sta97]. 
Here and and below, g is a variable, and not a complex number. The fundamental theorem of Wilf-Zeilberger 
states that a multi-dimensional finite sum of a (proper) g-hyper-geometric term is always q-holonomic; see 
[WZ92, ZeiQO; PWZ96]. Given this result, one can easily generate g-holonomic sequences. We learnt about 
this astonishing result from Doron Zeilberger in 2002. Putting this together with the fact that many state-sum 
invariants in Quantum Topology are multi-dimensional sums of the above shape, it follows that Quantum 
Topology provides us with a plethora of g-holonomic sequences of natural origin; [GL05]. For example, 
the sequence of Jones polynomials of a knot and its parallels (technically, the colored Jones function) is 
g-holonomic. Moreover, the corresponding minimal recursion relation can be chosen canonically and is 
conjecturally related to geometric invariants of the knot; see [Gar04]. Recently, the author focused on the 
degree of a g-holonomic sequence, and in the case of the Jones polynomial it is also conjecturally related to 
topological invariants of knot complement; see [Garll]. Since little is known about the Jones polynomial of 
a knot and its parallels, one might expect no regularity on its sequence of degree. The contrary is true, and 
in fact is a property of g-holonomic sequences and the focus of this paper. Our results were announced in 
[Garll] and [Gar], where numerous examples of geometric/topological origin were discussed. 

1.2. The degree and the leading term of a q-holonomic sequence. Our main theorem concerns the 
degree and the leading term of a g-holonomic sequence. To phrase it, we need to recall what is a g-holonomic 
sequence, and what is a quasi-polynomial. A sequence (/n(g)) of rational functions is q-holonomic if it 
satisfies a recursion relation of the form 

(1) a^q", q)fn+d{q) + ■■■ + ao(g", q)fn{q) = 

for all n where aj{u, v) G Q[u, v] for j = 0, . . . , d and ad{u, v) ^ 0. 



Given a polynomial with rational coefficients 

M 

L 



/(g)=Vc,g' 



that satisfies CmCM 7^ 0, its degree (also known as the order in different contexts) 6{f(q)) and its leading 
term lt[f[q)) is given by m and Cm respectively. In other words, the degree and the leading term of f{q) is 
the order and the starting coefficient in the Taylor series expansion of /(g) at g = 0. The degree and leading 
term can be uniquely extended to the field Q(g) of all rational functions, the ring Q[[g]] of formal power 
series in g, the field Q((g)) of all Laurent series in g and finally to the algebraically closed field 

(2) K^Q{U}}^^r=Miq'^l) 

of all Puiseux series in g with algebraic coefficients (see [Wal78]). The field K is required in Theorems 1.2 
and 1.3 below. 

Recall that a quasi-polynomial pin) is a function 

d 
p : N — > N, pin) = ^ Cj{n)n^ 

3=0 

for some d G N where Cj{n) is a periodic fuction with integral period for j = 1,. . . ,d; [Sta97, BR07]. If 
Cj ~ for j > 2, then we will say that p is a quadratic quasi-polynomial. We will say that a function is 
eventually equal to a quasi-polynomial if they agree for all but finitely many values. 
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Theorem 1.1. (a) The degree of a g-holonomic sequence is eventually a quadratic quasi-polynomial, 
(b) The leading term of a g-holonomic sequence eventually satisfies a linear recursion with constant coeffi- 
cients. 

Theorem f .1 follows from Theorems 1.2 and 1.4 below. 

Remark 1.1. If fn{q) G '^[q^^] is g-holonomic with degree (5„ and leading term lt„, it follows that fn{q) — 
lt„g''" is also g-holonomic. Thus, Theorem 1.1 applies to each one of the terms of fn{q)- 

Remark 1.2. L. Di Vizio brought to our attention that results similar to part (a) of Theorem 1.1 appear 
in [BB92, Thm.4.1] and also in [DV08, Sec. 1.1, See. 1.3]. However, the statement and proof of Theorem 1.1 
appear to be new. 

Our next corollary to Theorem 1.1 characterizes which sequences of monomials are g-holonomic sequences. 
In a sense, such sequences are building blocks of all g-holonomic sequences. 

Corollary 1.3. If (a„) and {b„) are sequences of integers (with a„ ^ for all n), then (a„g''") is g-holonomic 
if and only if &„ is holonomic and a„ is a quadratic quasi-polynomial for all but finitely many n. 

1.3. The Newton polygon of a linear q-difTerence equation. As is common, we can write a linear 
q-difference equation (1) in operator form by introducing two operators L and Af which act on a sequence 
ifniq)) by 

(LfUq) = /„+!((?), {MfUq) = g«/„(g). 

It is easy to see that the operators M and L g-commute 

(3) LM = qML 
and generate the so-called q- Weyl algebra 

(4) V = U^iA'((A/i/'-))(i)/(LA/i/'' - q^Z-'ML). 

We will call an element of 2? a linear q-difference operator. The general element of V is of the form 

d 

(5) P = Y,a^{M,q)L' 

where a,{M, q) e K{{M^/'')) for i = 0, . . . , d and ad{M, q) ^ 0. The equation 

P/ = 

for a A"-valued sequence fn{q) is exactly the recursion relation (1). The Newton polygon N{P) of an element 
P € I? is defined to be the lower convex hull of the points (i, duiO'i)) where 5m denotes the smallest degree 
with respect to M . Note that usually the Newton polygon N'{P) of a 2-variable polynomial is defined to 
be the convex hull of the exponents of its monomials. Since we are working locally at g = 0, we view N' [P) 
by placing our eye at — oo in the vertical axis and looking up. The resulting object is the lower convex hull 
N{P) defined above. The Newton polygon of a linear g-difference equation was also studied in the recent 
Ph.D. thesis by P. Horn; see [Hor09, Chpt.2]. 

The Newton polygon N{P) of P is a finite union of intervals with rational end points and two vertical 
rays. Each interval with end-points («, di) and (j, dj) for i < j has a slope s = (dj — di)/ {j — i) and a length 
(or multiplicity) I = j — i. The multiset of slopes s{N{P)) of N{P) is the set of slopes of s N{P) each with 
multiplicity Is- An example of a Newton polygon of a linear q-difference operator of degree 7 with three 
slopes —1, 0, 1/2 of multiplicity 2, 1, 4 respectively is shown here: 
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The Newton polygon is a convenient way to organize solutions to a linear q-difference equation. The reader 
may compare this with the Newton polygon of a linear differential operator attributed to Malgrangc and 
Ramis; see [vdPSOS, Sec. 3. 3] and references therein. In analogy with the theory of linear differential operators, 
we will say that P is a regular- singular g-difference operator if its Newton polygon consists of a single 
horizontal segment. In other words, after a minor change of variables, with the notation of (5) this means 
that a,{M, q) £ K[[M^/'']] for all i and ao(0, 17)0^(0, q) ^ 0. 

1.4. WKB sums. Since Equation (1) involves two independent variables q and g", let us set q" = u and 
consider the (/-difference equation 

(6) ad{u, q)f{uq'^, q) ^ h ao(u, q)f{u, (?) = 

Note that if f{u,q) solves (6), and if fn{q) := /('/"i'?) makes sense, then /„(g) solves (1) and vice-versa. It 
turns out that the general solution to (6) is a WKB sum. The latter are generalizations of the better known 
generalized power sums discussed in Section 3.1 below. 

Definition 1.4. (a) A formal WKB series is an expression of the form 

(7) Mu,q) = q^"\{qrA{n,u,q) 
where r = (7, X{q),A) and 

• the exponent 7 is a rational number, 

M 00 

A{n, u, g) = 5m (|>^,k{Q>''^'''n' 

i=0 fc=0 

is a polynomial in n with coefficients in Q((9^^'^))[[u^''^]] for some r G N. 

• A(9)eQ((gi/'-))and0o,o(g) = l. 

• there exists c G Q so that 5{(f)i,k{q)) > ck for all i and k 

(b) A formal WKB sum is a finite K-\mea.r combination of formal WKB series. 
Observe that if fr{u,q) is a formal WKB series, then its evaluation 

(8) frA'l)^fr{q",q)^K 

is a well-defined i^T-valued sequence for n > —re. 

Observe further if /t('u, q) is given by (7), the operators L and M act on friu, q) by: 

(9) {Lf^){u,q) ^ q-<'^^^+''>\{qr+'A{n + l,uq,q), iMf^){u,q) = q-'"\{qruA{n,u,q). 

Theorem 1.2. Every linear q-dijferenee equation has a basis of solutions of the form fT,ni<l)- Moreover, the 
multiset of exponents of the bases is the multiset of the negatives of the slopes of the Newton polygon. 

Recall that iir{{Af }} denotes the field of Puiseux series in a variable M with coefficients in K. Let (5m(/) 
denote the minimum exponent of M in a Puiseux series / G K{{M}}. 

Theorem 1.3. Every monic q-difference operator P £ T? of order d ean be factored as an ordered produet 

d 

(10) P^\{{L-a,{M,q)) 

4=1 

where ai{M,q) G K{{M}} and the multiset {SM{o-i)\i = 1, . . . ,ii} of slopes is the negative of the multiset of 
slopes of the Newton polygon of P. 

Remark 1.5. The WKB expansion of solution of linear ^-difference equations given in Theorem 1.2 appears 
to be new, and perhaps it is related to some recent work of Wittcn [Wit], who proposes a categorification of 
the colored Jones polynomial in an arbitrary 3-manifold. The curious reader may compare [Wit, Eqn.6.21] 
with our Theorem 1.2. We wish to thank T. Dimofte for pointing out the reference to us. 
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Remark 1.6. Although the proofs of Theorems 1.2 and 1.3 follow the well-studied case of linear differential 
operators in one variable x, we are unable to formulate an analogue of Theorems 1.2 and 1.3 in the differential 
operator case. Indeed, g is a variable which seems to be independent of the spacial variables x of a differential 
operator. The meaning of the variable q can be explained by Quantization^ or from Tropical Geometry (where 
it is usually denoted by t) or from the representation theory of Quantum Groups; see for example [Gar] and 
referencies therein, and also Section 4 below. 

Theorem 1.4. If fniq) is a finite K-linear combination of formal WKB series of the form (7), then for 
large n its degree is given by a quadratic quasi-polynomial and its leading term satisfies a linear recursion 
relation with constant coefficients. 

1.5. An example. In this section we discuss in detail an example to illustrate the introduced notions and 
the content of Theorems 1.1, 1.2 and 1.4. The example shows that the proof of Theorem 1.1 and the 
WKB expansion of Theorem 1.4 is algorithmic, despite the use of differential Galois theory. Consider the 
g-holonomic sequence fniq) G Z[g] whose first few terms are given by: 



/o -- 


= 1 


/l -- 


- 2-q' 


/2 -- 


= 3 + g - 2q^ - 2g^ + g^ 


/3 -- 


= 4 -f 2g + 2g" - 3g^ - 4g' - 4g'' - g V 2g° + 2q^° + 2g" - g^' 


h -- 


= 5 + 3g -f 4g^ -f 3g^ -f g'' - 4g" - 6q'^ - 8g^ - 8g* - 4g^ - 2g''' + 3g'' + 4g'' + 7g" + 5g" 




+4g^^ + q'' - 2g^« - 2g^« - 2g^° - 2q'' + q'" 


/5 = 


= 6 + 4g + 6g^ + 6g^ + 6q'^ + 2q - Zq - %q - \2q - ISg'' - 16g^" - llg" - 8g^^ + Sg^"* 




+ 12g^^ + 14g« + 16g^^ + 12g^« + 10g^« + 4g^° - q^^ - 4g^^ - 7q'' - 8g^^ - 8g^^ - 5g^« - 4g^^ 




-q'^ + 2q^° + 2q'^ -f 2q^^ + 2g^^ -f 2g^* - g« 


/6 = 


= 7 + 5g + 8g^ -h 9g^ -h llg"* + Qq' + 7g'' - 2g' - 7g** - ISg** - 22g" - 28g" - 30g'' - 26g" 




-22g" - llgis - 2q"' + Uq^' + 21q^^ + 33g"^ + 34g2" + 36g^^ + 30g"^ + 25g^^ + llg^^ + 3g^^ 




-iq'" - nq" - 22g^« - 24g^^ - 24g^" - 20g«^ - lAq'^ - IQq'' - q'^ + 2q'' + 7q''' + 8g«" 




+llg^** + 9q"' + 8g" + 5g« + 4q^^ + q« - 2g« - 2g« - 2g^^ - 2g^« - 2g^« - 2g'^° + g" 



The general term of the sequence / is given by: 
(11) W'Z)=E7-f^ 

^ (q)n-k(q)k 

where the q-factorial is defined by 

(g)n=n(i-9'). (?)0=1- 

fe = l 

The above expression implies via the Wilf-Zeilberger theorem that / is q-holonomic; see [WZ92]. The qzeil.m 
implementation of the Wilf-Zeilberger proof developed by [PR97, PR] gives that / is annihilated by the following 
operator: 



Pf = 



(-1 + Mq^){-1 + Mq + M^q^)L^ + (-2 + Mq + Mq^ + 2M^g^ + M^q^ + 2M^q* - 2M^g* - 
2M=*g'' - M^q^ - M^q'' - M\'' + M^'q' + M\'' + Afq')L + 1 - Mq^ - Af^g^ 
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The 2-diinensional Newton polytope M2{Pf) and the Newton polygon N{Pf) are given by: 

MA A 



mPf)^ 




Pf is a second order regular-singular g-difference operator. Its Newton polygon N{Pf) has only one slope with 
multiplicity 2. The edge polynomial of the 0-slope is cyclotomic (L — 1)^ with two equal roots (i.e., eigenvalues) 1 
and 1. The degree 5„ and the leading term lt„ of fn{q) are given hy 5„ — and lt„ = n + 1 for all n. 

The characteristic curve ch/ (discussed in Section 4) is reducible given by the zeros {L,M) £ (C*)^ of the 
polynomial 

(-1 + M + A/^)(-l + 2L - L^ + L^M - 3LAf + LAl'-^ + LM^) = 
The tropical . lib program of [Mar] computes the vertices of the tropical curve 7/ are: 

(-1,-2), (3,-2), (0,-3/2), (1,-3/2), (0,-1) 

The next figure is a drawing of the tropical curve 7/ and its multiplicities, where edges not labeled have multiplicity 
1. 




The Newton subdivision of the Newton polytope N2[Pf) is: 




To illustrate Theorem 1.1, observe that ,fn{q) is a sequence of Laurent polynomials. The minimum (resp., maxi- 
mum) degree 5„ (resp., (5„) of /n(g) with respect to q is given by: 



5n=0, 



5n — 



n{3n + l) 



The coefficient lt„ (resp., li„) of g " (resp., g " ) in /n(g) is given by: 

lt„=n-hl, lt„ = (-1)" 

lt„ and lt„ are holonomic sequences that satisfy linear recursion relations with constant coefficients. 
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To illustrate Theorem 1.4, observe that there is a single slope of length 2 with edge polynomial [L — 1)^ and 
eigenvalue 1 with multiplicity 2. This is a resonant case. Theorem 1.4 dictates that we we substitute the WKB ansatz 



+j)k 



k = 

in the recursion Pjf = 0, collect terms with respect to M = q and with respect to n and set them all equal to zero. 
It follows that the vector {tpk{q),4>k{q)) satisfies a sixth order linear recursion relation: 

^fe = T. ^ Yt'^'''4>-6+kq^*^{l + q)ip-5+k + q^^'' (! + <? + g^) ^-4+fc + g"^+'' [2q^ + 1q^ - q^ ^-3+fc 

/ 4 2fe-2 I o fe I 1 + fe I o 2+fe I -l+2fe\ , I / 2 I 2fe-l fe 1+fe , 2fc\ / 1 

- ( -g -q +2g+g^-|-2g^+g ^ j V)_2+fe + ( g +g -g-g^+g 1 V>_i+fc| 

'?^'= = ... _ fc^^2 {-g^"^''>-6+fc - g^"^''(l + q)(l>-5+k + g^"^*" (1 + g + g^) </)-4+fe + g^^"^*" (2g-^ + 2q* - g* j (/)-3+fc 
(-g« - g^'' + 2q2+'-- + g^+^- + 2g^+^- + q'+^'') 0_2+fe (g^ + q^" - g^+*^ - g^+'= + q'+"' 



,. , , , , , , i+fc 



+ 



g 

9^+^° (f + g^) i,-e,+k q'+\l + g) (1 + g'-) ^s+k q'^" (1 + g + g^) (l + g'') ^-4+^ 
(-l + g'^) (-i + g'=) (-1 + g*^) 

29-^+" (g^ + g^ - q" + q'+'' + q^+'') ^_3+fc 

(-i + g''-) 

^-2+k (2g2 + g3 ^ 2g^ - 2g« - 2q'' + 2q'+'' + 2q^+>' + g3+* + 2g'*+*) ^_2+, 



+ 



(-1 + gfc) 

g-^+'° (g + g^ - 2g^ - 29" - 9^+^^ + g^+'') ^_i+, 

(-1 + gfc) 



with an arbitrary initial condition {ipo, 4>o) G Q[[g]]^- Note that ipk the first equation above is a recursion relation for 
■i/ifc and the second equation is a recursion for 0^ that also involves ip^i for k' < k. This is a general feature of the 
WKB in the case of resonanse. It follows that our particular solution (11) of the g-difference equation Pff = has 
the form: 

00 

(12) /n(g) = J2(Mq) + nMq))q"'' e QIM 

fe=0 

where {ipic,4'k) are defined the above recursion with suitable initial conditions. Note that Equation (12) implies that 
the coefficient Cm,n = cim + nbm of g™ in fn{q) is a linear function of n for ?i > m. These values determine our initial 

conditions by: 

00 

00 (g) = ^ flmg'" 

m = 

and 

^o(g) = ^ fcmg™' 

In fact, a direct computation shows that the first 31 values of am for m = 0, ... ,30 are given by: 

1, -1, -4, -9, -19, -33, -59, -93, -150, -226, -342, -494, -721, -1011, -1425, -1960, -2695, -3633, -4903, 

-6506, -8633, -11312, -14796, -19157, -24773, -31744, -40608, -51578, -65372, -82341, -103522 
and the first 31 values of bm for 771 = 0, . . . , 30 are given by: 

1, 1, 2, 3, 5, 7, 11, 15, 22, 30, 42, 56, 77, 101, 135, 176, 231, 297, 385, 490, 627, 792, 1002, 1255, 1575, 

1958,2436,3010,3718, 4565,5604 
The reader may recognize (and also confirm by an explicit computation) that 



1 °° 1 

g;g 00 -^-^, 1 -g 
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An extra-credit problem is to give an explicit formula for the power series (f>o{<l)- Note finally that Equation (11) 
implies that the specialization of /„ (q) at g = 1 is given by: 

fn{l) = 1 

for all n, much like the case of the colored Jones polynomial of a knot, [GL05]. 

1.6. Plan of the paper. Theorem 1.1 follows from Theorems 1.2 and 1.4. 

Theorem 1.2 follows by two reductions using a g-analogue of Hensel's lemma, analogous to the case of linear 
differential operators. The first reduction factors an operator with arbitrary Newton polygon into a product of 
operators with a Newton polygon with a single slope. After a change of variables, we can assume that these operators 
are regular-singular. The second reduction factors a regular-singular q-difference operator into a product of first 
order regular-singular operators with eigenvalues of possibly equal constant term. Thus, we may assume that the 
g-difference operator is regular-singular with eigenvalues of equal constant term. In that case, we can construct a 
basis of formal WKB solutions of the required type. The above proof also factorizes a linear g-difference operator 
into a product of first order g-difference operators, giving a proof of Theorem 1.3. 

Theorem 1.4 follows from the Lech-Mahler-Skolem theorem on the zeros of generalized exponential sums. 

Finally, in Section 4, which is independent of the results of the paper, we discuss several invariants of g-holonomic 
sequences, and compute them all in some simple examples. 

2. Proof of Theorem 1.2 

2.1. Reduction to the case of a single slope. Theorem 1.2 follows ideas similar to the case of linear differential 
operators, presented for example in [vdPSOS, Sec. 3]. Rather than develop the whole theory, we will highlight the 
differences between the differential case and the g-difference case. 

In the case of linear differential operators, the basic operators x and xdx satisfy the inhomogeneous commutation 
relation 

(13) [xdoojx] = X 

The right hand side of the above commutation relation leads to consider Newton polygons are convex hulls of translates 
of the second quadrant R_ x R^. For examples of such polygons, see [vdPS03, Sec. 3. 3]. As a result, the slopes of 
the Newton polygons of linear differential operators are always non-negative rational numbers. On the other hand, 
in the g-difference case, the g-commutation relation (3) preserves the L and M degree of a monomial in M, L. 

In the case of linear differential operators of a single variable x, the WKB solutions involve power series in x^''^ 
and polynomials in logx. In the case of linear g-difference operators, the WKB solutions involve power series in g" 
and polynomials in n. In addition, the convergence conditions are easier to deal with when g = 0, and much harder 
when g = 1. Convergence when g = 1 involves regularizing factorially divergent formal power series. 

Consider a g-difference operator P and its Newton polygon N{P). In this section, we will allow the coefficients 
ai{M, g) of P from Equation (5) to be arbitrary elements of the field K{{M}}. Recall that the Newton polygon N{P) 
consists of finitely many segments and two rays. We will say that P is slim if all its monomials lie in the bounded 
segments of its Newton polygon. We will say that Pi > P2 if the boundary of N{Pi) lies in the interior of P2 or in 
the interior of the two rays of P2. We will say that an operator P from Equation (5) is momc if ad{M, g) = 1. Recall 
the Minkowski sum 

Ni+N2 = {a + b\ae Ni, b e N2} 

of two polytopes in the plane. Here and below, we will always consider the Minkowski sum of two convex polytopes 
that arise from linear g-difference operators, i.e., polytopes that consist of two vertical rays and some bounded 
segments. 

The next lemma is a g- version of Hensel's lifting lemma. 

Lemma 2.1. If P is monic and slim, for every partition of its set of slopes in two different sets Si and S2 there exist 
unique slim and monic operators Pi and P2 with slopes Si and S2 respectively, and an operator R{P) so that 

P = PiP2 + P(P), R{P)>P. 

Proof. This follows as in the differential case. The only difference is that monomials M'^L^ and M" L!' commute up 
to a power of g. The next example will explain the proof. Suppose P — aM + 6i + cM^L + M'^L^ is a monic slim 
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operator with Newton polygon N with two slopes, where a,b,c £ K. Now N = Nj + N2 shown as follows: 







+ 




We want to find unique sq, yo, J/i, J/2 £ K so that 



aM + bL + cKPL + M'-^L^ = {xoM + L){yo + yiML + y-zM^L^). 



Multiplying out, we obtain that 
aM + bL + cA'fL + N'fL^ = , 
Thus, we obtain the system of equations 



aM + bL + cK'fL + A/^L^ = xoyoM + yoL + yiLML + R{P), R{P) = xoyihf L + xqM^L'^ + LKfL^ > P. 



xoyo = a, 2/0 = 6, qyi = c. 

The above system can be solved uniquely in K whether q is present or not. 



D 



Proposition 2.2. Suppose P is monic and N(P) = Ni + N2 where Ni and N2 have no slope in common. Then 
there exist unique momc operators Pi and Pz with Newton polygons Ni and N2 such that P = P1P2. In addition, we 
have: 

V/VP = V/VPi ® V/VP2. 



Proof, ft follows verbatim from the proof of Theorem 3.48 in [vdPSOS] using Lemma 2.1. 



D 



2.2. Reduction to the case of a single eigenvalue. Suppose that P is a linear g-difference operator with a single 

2 I 
slope. After a change of variables fn{q) = ?" ^ '5'?i(?), we can assume that the slope is zero, i.e. that P is regular 

singular. In other words, if P is given by (6), then 0^(0, g)ao(0, q) 7^ 0. Then, the slim part S{P) of P is given by 

d 

S(P) = ^a,(0,g)i' 

i=0 

where P — S{P) > P. We may think of S{P) as a polynomial in a variable L with coefBcients in the field K. Suppose 
that S{P) is monic and we can factor S{P) = AB as the product of two relatively prime monic polynomials. 

Proposition 2.3. With the above assumptions there exist unique monic regular-singular operators Pi and P2 such 
that P = P1P2 and S{Pi) = A and S{P2) = B. Moreover, 

V/VP = V/VPi ® V/VP2. 
Proof. It follows verbatim from the proof of Proposition 3.50 in [vdPS03] using Lemma 2.1. □ 

2.3. First order linear g-difference equation. In this section we will prove Theorem 1.2 for the first order linear 
g-difference equation 

fn+i{q) -b{q",q)fn{q) = 
where b{M,q) € K{{M)). If b{M,q) = M^c{hd,q) where c{M,q) G K[[M]] and c(0, g) 7^ the substitution /„(g) = 



c(0, g)"(7n(g) gives the equation 



(14) 



5n+i(g) -a{q",q)9n{q) = 



where a{M,q) = q'^^^^^ G K[[M]] satisfies a(0, g) = 1. A solution to Equation (14) is 

fn[q) = a{q,q)a{q'^,q)...a(q"''^,q) 

nr=i«(9'.9) 



Now 



a(g",g)a(g"+i,g)... 
a(M,g) = l-^M'=afc(g) 
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where ak{q) £ K, implies that 



oo / oo 



a 



{M,q) 



TTT = 1 + E E^'«^W (^^K[[u]] 



Consequently, 

Let us now match this solution with a WKB sum. Consider 

oo 
fc=0 

where <j>o = ^ and substitute into the equation 

f{uq, q) - a{u, q)f{u, q) = 
we obtain an equation in the ring /('[[ii]]. The vanishing of the coefficient of u'' gives 



k 
i 

Ok 



1 * 



1 .=1 



Thus the WKB solution matches exactly the solution (15) for the first-order regular singular equation (14). 

2.4. Proof of Theorem 1.2. Propositions 2.2 and 2.3 and a possible replacement of q by q^^^ reduce Theorem 1.2 
to the case of a monic regular-singular differential operator of the form 

d 

(16) P^l[{L-a,{A4,q)) 

!=1 

where ai{M, q) € /('[[A/]] for all i and a(0, q) — 1. It suffices to find a basis of solutions of the equation Pf = of the 
form 

oo 

/(">'?) = '^4>k{n,q)u'' 

k = 

where (f)k(n, q) £ K[n] are polynomials of degree at most d. We will use induction on d. When d = 1, this is proven 



in Section 2.3. Let Q = Yli^ii^ ~ "-ii^, q))- If we set gn{q) = fn+i{q) — ai{q", q)fn{q) then it Pf = implies that 
05 = 

(17) 



It 

Qg = 0. Conversely, the system Pf = is equivalent to the system of equations 

f{uq,q) -ai{u,q)f{u,q) = g{u,q) 
(u,g) = 

Now, by induction we have 

oo 

9{u,q) = 'pipkin, q)u' 

fc=0 

where ipkin, q) € K[n] are polynomials of degree at most d — 1. Set f{u, q) — X^feLo 4>k{n, q)u^ where 4)k{n, q) £ K[n\ 
are polynomials of degree at most d in the first Equation (17). Both sides are elements of /("[[ujjfn]. The vanishing of 
each coefficient of u ri computes the (j)k in terms of the ijik by a rank 1 linear system of equations. It follows that we 
can find a d-dimensional linear space of WKB series solutions of Pf — 0. This concludes the proof of Theorem 1.2. 

2.5. The regular-singular non-resonant case. For completeness, we give a short proof of Theorem 1.2 in the case 
of a non-degenerate regular- singular operator. Recall from Section 1.3 that a regular- singular g-difference operator 
has the form 

d 

(18) P = ^a.(M,g)r 

i=0 

where ai{M,q) £ A'[[A/^''']] for all i and ao(0,q)ad{0,q) / 0. The characteristic polynomial x{^tM,q) of of a 
regular-singular operator P is given by 

d 

(19) x(a;,M,g)-^a,(Af,g)a;' 
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The eigenvalues of P are the roots of the equation xi^t ^^i q) ~ 0- Since K is algebraically closed, Puiseux's theorem 
implies that the eigenvalues of P are elements of the field -R'liq}}; see [Wal78]. In other words, a regular-singular 
operator P has d eigenvalues Xi{M,q) G K{{M^''^ )) for some r' . After we replace r by r', we assume r' — r below. 
We say that a regular singular operator P is non- degenerate if its eigenvalues satisfy the non-resonanse condition 
that Aj(0, q)/Ai(0, q) / q^ for aU i / j and aU k. 

Suppose that P is a non-degenerate regular-singular operator as above. Consider Equation (6) 

(20) ad{u, q)f{uq'^, q) -\ h ao{u, q)f{u, g) = 

We will show that the general solution to (20) is a WKB sum. Without loss of generality we assume r — 1. We will 
construct a basis of solutions of the Equation (20) of the form 

oo 

(21) /(«,?) = A(g)"^0fe(g)n'= 

fc=0 

that satisfy (^o(g) = 1- Here n = log it/ log g. The point is that although f{u, q) ^ K^u\\, the ratio f{uq-', g)/A(g)" £ 
isr[[u]] for every j. This follows from Equation (9). Substitute (21) into (6) and divide by A(g)". The result is an 
equation in iC'[[?i]]. The vanishing of the coefficient of u^ gives 

d 

^a,(0,g)A(g)^=0 
i.e., X — A(g) is a root of the equation x{^^ 0, g) = 0. The vanishing of the coefficient of u^ for fc > is the condition 

(22) 0fc(g)co,fe(g) + '^^cj)k-i{q)ci^k-i{q) = 

i = l 

where 

Ci,j{q) = (d/rfM)'x(a;,^, ?)!(:,, i,)=(5JA(9),o)- 
Using the initial condition 00 = 1, it follows inductively that 

fe 
i/'fe(g) ■■= <t)k{q)Y{co,j{q) 
j=i 

is a homogeneous polynomial of total degree k in the variables Cij{q), where the degree of Cij{q) is i + j. For example, 
we have: 

Tpi = -ci,o 

■02 = Ci,oCi,l — Co,lC2,0 

■03 = — Ci,oCl,lCl,2 + Co,lCl,2C2,0 + Co,2Cl,oC2,l — Co,lCo,2C3,0 

■04 = Ci,oCl,lCl,2Cl,3 — Co,lCl,2Cl,3C2,0 — Co,2Cl,oCl,3C2,l — Co,3Cl,oCl,lC2,2 + Co,lCo,3C2,OC2,2 

+Co,lCo,2Cl,3C3,o + Co,2Co,3Cl,oC3,l — Co,lCo,2Co,3C4,0 
05 = — Cl,oCi,iCi,2Cl,3Ci,4 -|- Co,lCl,2Cl,3Cl,4C2,0 + Co,2Cl,oCi,3Ci,4C2,l -|- Co,3Cl,oCl,lCi,4C2,2 

— Co,lCo,3Cl,4C2,OC2,2 + Co,4Cl,oCl,lCl,2C2,3 — Co,lCo,4Cl,2C2,OC2,3 — Co,2Co,4Cl,oC2,lC2,3 

— C0,lC0,2Cl,3Cl,4C3,() + CO,lCO,2Co,4C2,3C3,0 — Co,2Co,3Cl,OCl,4C3,l — Co,3CO,4Cl,oCl,lC3,2 
+C0,lCo,3C0,4C2,0C3,2 -f C0,lC0,2C0,3Cl,4C4,0 + C0,2C0,3C0,4Cl,0C4,l — C0,lC0,2C0,3C0,4C5,0 

Note that co.fe = x{l''Ml)y^y'l) 7^ ^^r ^■'1 ^ due to the non-degeneracy of (20). It remains to show that there 
exists c such that 5(0fc(g)) > ck for all k. To prove this, observe that 

co,fc(g) = X^a,(0,g)(g'=A(g))^ 

3=0 

satisfies that 

S{co,k) = mm{S{aj{0, g)A(g)^) + fcj} 

0<j<d 

is attained at least twice. Since ao(0, g) ^ 0, it follows that for large fc, the above minimum is attained at j = 0. In 
other words, we have 

S{co,k) = C-- 5{ao{0,q)). 
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Thus, 

fc 

S{Y[coAq))>ck + c 

for all k large enough. On the other hand, S{ci,j{q)) > for all i,j. The above formulas imply that 5{(f)k{q)) > c'k 
for all k. Thus, for every eigenvalue A(0, g) of P there is a unique WKB solution (21) of Equation (20). Since there 
are d distinct eigenvalues, the constructed solutions form a basis for the vector space of solutions of (20). This gives 
a proof of Theorem 1.2 in the non-degenerate regular-singular case. 

Remark 2.4. In the case of regular-singular linear differential equations, expressions of the form (22) appear. 
However, one needs to estimate the numerator of those expressions, as well as the denominator. See for example, 
[GG06] and references therein. This explains why the WKB theory of g-difference equations at q = is much simpler 
than the corresponding theory at g = 1. 

3. Proof of Theorem 1.4 

3.1. Generalized po^wer sums. An important special case of Theorem 1.1 is the case of a linear recursion with 
constant coefficients. In this rather trivial case, for every n, fn{q) is a constant function of q, so the degree is easy to 
compute. 

Generalized power sums play a key role to the SML theorem. For a detailed discussion, see [vdP89] and also 
[EvdPSW03]. Recall that a generalized power sum an for n = 0, 1, 2, . . . is an expression of the form 

m 

(23) a„ = ^yli(n)ar 

with roots cti, 1 < i < m, distinct nonzero quantities, and coefficients Ai{n) polynomials respectively of degree 
n{i) — 1 for positive integers n(i), 1 < i <m. The generalized power sum a„ is said to have order 

1=1 
and satisfies a linear recursion with constant coefficients of the form 

In + d = Slffln + d-l + ' ' ' + Sda^ 

where 

m 

s{x) = \{{1 - a,a;)"('' = 1 _ Six - . . . s^x" . 

It is well known that a sequence satisfies a linear recursion with constant coefficients if and only if it is a generalized 
power sum. The LMS theorem concerns the zeros of a generalized power sum. 

Theorem 3.1. [Sko35, Mah35, Lec53] The zero set of a generalized power sum is the union of a finite set and a 
finite set of arithmetic progressions. 

3.2. Proof of theorem 1.4. Consider a formal WKB series 

2 I r — ^ j + nk 

/r,n(g) = q" ^ ''"^" 2_^ Cr,j,k{n)q - 
j,k>0 

where rj = 5{X{q)), /i — lt(A(q)) G Q and CTj,k{n) are generalized power sums with rth roots of 1. Fix a sequence 
fn{q) given by a finite Q((q^'''"))-linear combination of formal WKB series fr-{u,q). It follows that 

(24) /„(q) = E^^ (?)/-.. -(<?) = E?"'^'^""/^" E c.,,,fc(n)q"^ 

iei iei j,k>o 

where 7 is a finite set, Ci(q) £ Q((q^ )) are non-zero and (r^, Ai) are pairwise distinct. Consider the subset I' of / 
where the minimum 

min{7,} 

is achieved. 

Case 1. If /' consists of a single element io and Cio,jo.'=o('^) i^ ^'^^ identically zero, then let us concentrate on the 
following part of f„ (q) : 

„{q)^q ''" ''<> fJ.igCigjg^koWq - +■■■ 
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The LMS theorem concludes that there is a finite set of fuU arithmetic progressions that cover the set of natural 
numbers with finite complement, such that the restriction of the degree of /n(?) to each full arithmetic progression 
(with the exception of a finite set) is given by 'n?'yi„ +r)ign+ if>+^'^° and the leading term is given by At"„Cip_jp.fc,-,(n). 
When Ci(jj(,,fcQ(n) vanishes consider the next term in the series (24) and repeat the above argument. After finitely 
many repetitions the process will stop since fTi,n{q) are Q((g^''^))-linearly independent. 

Case 2. If/' consists of more than one element, without loss of generality, assume that Cififi{n) are not identically 
zero for i £ I' . Consider the subset I" of /' where the maximum 

maxliji} 

is achieved. 

Case 2.1. /" consists of a single element iq. Then, it follows that 

/n(g) =?"'"'« +'"o'VroC,o,0.0('^) + ... 

and the proof proceeds as in Case 1. 

Case 2.2. J" consists of more than a single element io- In that case, it follows that 

ie/" 
The leading term of /n(g) is a generalized power sum and the proof proceeds as in Case 1. 

4. Invariants of q-holonomic sequences 

4.1. Synopsis of invariants. In this section, which is independent of the results of our paper, we summarize various 
invariants of a a g-holonomic sequence. Some of these invariants were announced in [Gar04, Garll, Gar]. In this 
section / denotes a g-holonomic sequence /„(?) £ Q(g) that satisfies Equation (1) where ai{A'I, q) £ Q(Af, q). Here is 
a summary of invariants of /. 

(a) The annihilator polynomial Pf 

(b) The 2 and 3-dimensional Newton polytopes. 

(c) The characteristic curve ch/. 

(d) The tropical curve Tf . 

(e) The degree and leading term of /. 

As was explained in [Gar], these invariants fit well together and read information of the g-holonomic sequence /. For 
completeness, we will present these ideas here, too. 

4.2. The annihilating polynomial of a g-holonomic sequence. In this section we discuss the annihilating 
polynomial of a g-holonomic sequence /. Informally, it is the minimal order recursion relation (homogeneous or not) 
that / satisfies. This was studied in [Gar04] and [Gar]. We will call an operator P given by (5) with r = 1 reduced if 
ai{M,q) £ Z[M, g] are coprime. Let 

V = Q(A/, q){L)/{LM - qML) 

denote the localized g-Weyl algebra. 

First, we define the homogeneous annihilator Pf of /. Consider the O'-module (nonzero, since / is g-holonomic) 

(25) Mf ^{P eV' \Pf = 0} 

and its unique monic generator P' G D'. After multiplying by the common denominator of the coefficients of P' with 
respect to the powers of L, we obtain a unique reduced generator Pj. 

Next, we define the non-homogeneous annihilator {Pf,bf) of /. Consider the O'-module 

(26) Mf = {{P, b)eV' X q{M, g) \Pf = b} 

Mf has a unique generator of the form {P',e) where P' £ V' is monic (with respect to L) and e £ {0, 1}. After 
multiplying by the common denominator of the coefficients of P' with respect to the powers of L, we obtain a unique 
generator of the form (Pp , 6/) where Pf' is reduced and bf £ Z[M, gj. 

Definition 4.1. We call Pf and [Pf, bf) the homogeneous and the non-homogeneous annihilator of the g-holonomic 
sequence /, respectively. 

In [Gar, Lem.1.4] we show how Pf determines {Pf,bf) and vice-versa. 
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4.3. The characteristic variety of a q-holonomic sequence. The characteristic curve of a q-holonomic sequence 
/ is the zero set 

(27) ch^ = {(L,M) G (C*)^ |P^(i,M,l) =0} 

In case / is the colored Jones polynomial of a knot K, the AJ Conjecture of [Garll] identifies the characteristic curve 
with a geometric knot invariant, namely the moduli space of SL(2,C) representations of the knot complement. 

4.4. The Newton polytope of a g-holonomic sequence. In this section we consider the 3-dimensional and the 
2-dimensional Newton polytope of a g-holonomic sequence /. Let us write the annihilating polynomial P/ of / in 
terms of monomials: 

(28) Pf= Y. '^^,3,k'l'M'L\ 

[i,j,k)£A 

where a,,j,fe G Q and the sum is over a finite subset ^ of N which depends on /. 

Definition 4.2. (a) The 3-dimensional Newton polytope Af-^{Pf) of / is the convex hull A/'3(P/) of A. 

(b) The 2-dimensional Newton polytope 7V2(P/) of / is in R^ is the projection of A/'3(P/) on R^ under the map 

{x,y,z) !-)■ {x,y). 

The width of A/3 (P/) with respect to L is the degree of a minimal order recursion relation for /, and it is a measure 
of the complexity of /. Likewise, the width of A/'3(P/) with respect to M is a measure of the complexity of the 
coefficients of a recursion relation of /. Knowing (or guessing) Afs{Pf) given some terms in / has applications in 
Quantum Topology. 

Generically one expects that the Newton polygon of the polynomial Pf(x,y, 1) coincides with 7V2(P/). 

Lemma 4.3. The lower convex hull of A/'2(P/) coincides with the Newton polygon N{Pf) of Pj from Section 1.3. 

4.5. The tropical curve of a g-holonomic sequence. In this section, we assign a tropical curve 7/ to a g- 
holonomic sequence /, following [Gar]. For a leisure introduction to tropical geometry, see [RGST05, SS09, Stu02, Gar] 
and references therein. The main idea is to think of an operator in M, L with coefficients in Q{q) as a g-parameter 
family of polynomials in two commuting variables. More precisely, the coefficients of the annihilating polynomial Pf 
given by (28) define a function 

(29) P:R^ — ^R, F(x,y)= min {ix + jy + k} 

{i,i,k)eA 

P is a piecewise linear convex function. The locus of the points in R'^ where P is not differentiable is the tropical 
curve Tf of /. By definition, 7/ is the locus of points (a;, y) G R^ where the minimum in (29) is achieved at least twice. 
It is well-known that 7/ consists of a finite collection of line segments with rational vertices, and a finite collection of 
rays with rational slopes, together with a set of multiplicities that satisfy a balancing condition. Abstractly, a tropical 
curve is a balanced rational graph, and vice- versa; see [RGST05, Thm.3.6]. 

There is a duality between Newton polytopes and tropical curves. Indeed, the projection of the lower hull of 
Mi{Pf) gives a Newton subdivision oi J\f2(Pf). The tropical curve 7} is dual to the Newton subdivision of A/'2(P/). 
For drawings of tropical curves of g-holonomic sequences of geometric origin, see [Gar]. 

4.6. A tropical equation for the degree of a g-holonomic sequence. In this section we explain the relation 
between the degree 5{n) — 5(/„(g)) of a g-holonomic sequence / and its tropical curve 7/, following [Gar]. Since / 
is annihilated by Pf (given by (28)), it follows that for all natural numbers n we have: 

(i,j,k)eA 

Divide both sides by fn{q) and look at the degree with respect to q. It follows that for all ?i, 5{n) satisfies the 
following tropical equation: the minimum 

(30) min {5{n + i) - 5{n) + jn + fc} 

(i,j.k)&A 

is achieved at least twice. Although (30) may have non quasi-polynomial solutions. Theorem 1.1 suggests to look for 
quadratic quasi-polynomial solutions of (30). It turns out that those solutions 5{n) are determined by the tropical 
curve 7/ as follows. For n in an arithmetic progression minus finitely many values, we have: 

(31) S{n) = —n -fcin-fco 



E 
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Proposition 4.4. (ci,C2) satisfy the system of equations 

(32) C2i+j = C2i' + j' 

(33) -C2i + cii + k — -C2i' + cii' + fc' 
for two distinct points {i,j, k) and (i , j , k ) of A- 

Proof. For n in an arithmetic progression minus finitely many values we have: the minimum 

(34) min {n(c2i + j) + -C2i^ + cii + k} 

(i,j,k)eA 2 

is achieved at least twice. It follows that there are two distinct points {i,j, k) and {i' ,j' , k') of A such that Equations 
(32) and (33) hold. Equation (32) involves only (i, j) and {i' ,j') and says that C2 is the negative of a slope of N{Pf). 
Likewise Equation (32) involves only {i,k) and {i',k'). D 

Remark 4.5. It would be interesting to study the solutions to the tropical equation (30) independent of differential 
Galois theory. 

Remark 4.6. Each Equation (32) and (33) describes a tropical curve in R . However, the system of Equations (32) 
and (33) is not the intersection of the two tropical curves since the two planes R'^ are different projections of R'*. 

4.7. Future directions, g-holonomic sequences of several variables fni,...,nri<l) £ '^[q ] bIso appear in Quantum 
Topology. For example the colored Jones function of a 2-component link, of the s la-colored Jones polynomial of a 
knot is a g-holonomic sequence of two variables. Suitably formulated, Theorem 1.1 should extend to g-holonomic 
sequences of many variables. This will be explained in a future publication. 
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